Measurement of persistence in 1-D diffusion 
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Using a novel NMR scheme we observed persistence in 1-D 
gas diffusion. Analytical approximations and numerical sim- 
ulations have indicated that for an initially random array of 
spins undergoing diffusion, the probability p{t) that the av- 
erage spin magnetization in a given region has not changed 
sign (i.e., "persists") up to time t follows a power law t~^ , 
where 9 depends on the dimensionality of the system. Using 
laser-polarized '^■^^Xe gas, we prepared an initial "quasiran- 
dom" ID array of spin magnetization and then monitored the 
ensemble's evolution due to diffusion using real-time NMR 
imaging. Our measurements are consistent with analytical 
and numerical predictions of 6 « 0.12. 
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The dynamics of non-equilibrium systems is a field of 
great current interest, including such topics as phase 
ordering in binary alloys, uniaxial ferromagnets, and 
nematic liquid crystals, as well as coarsening of soap 
froth and diffusion of inhomogeneous fluids [e.g. 
The evolving spatio-temporal structures in these non- 
equilibrium systems depend crucially on the history of 
the system's evolution and are not completely charac- 
terized by simple measures such as two-time correlation 
functions. Therefore, an important problem in the study 
of non-equilibrium dynamics is the development of sim- 
ple and easily measurable quantities that give nontrivial 
information about the history of the system's evolution. 
The recently identified phenomenon of "persistence" may 
be such a quantity: it characterizes the statistics of first 
passage events in spatially extended non-equilibrium sys- 
tems P-pH . Persistence is being actively studied in sta- 
tistical physics; e.g., in the search for universal behavior 
in non-equilibrium critical dynamics . Practi- 

cally, persistence may be important in determining what 
fraction of a system has reached a threshold condition as 
a function of time; e.g., in certain chemical reactions or 
disinfectant procedures. 

Consider a non-equilibrium scalar field 0(x, t) fluctuat- 
ing in space and time according to some dynamics (e.g., 
a random array of interdiffusing spins). Persistence is 
the probability p{t) that at a fixed point in space the 
quantity [(/)(x, t) — (0(x, i))] has not changed sign up to 
time t. It has been found that this probability decays 
as a power law pit) ~ , where the persistence expo- 
nent 9 is generally nontrivial and has been suggested as a 



new universal dynamic critical exponent [plJl l|] . This ex- 
ponent depends both on the system dimensionality and 
the prevalent dynamics |^,^ , and is difficult to determine 
analytically due to the non-Markovian nature of the phe- 
nomena. Although 9 has been calculated - largely using 
numerical techniques ~ for such systems as simple diffu- 



sion §,|,|T|, the Ising model i,! 



15[, and the more gen- 
few measurements of 



eralized g-state Potts model |^,^ 
persistence have been performed (see Table In par- 
ticular, "breath figures" 2-D soap froth and 
twisted nematic liquid crystals are the only systems 
for which experimental results have been reported. Fur- 
ther experimental investigation is needed to probe the 
utility of persistence in studies of fundamental and ap- 
plied non-equilibrium dynamics. 

In this paper we present the first measurement of per- 
sistence in a system undergoing diffusion (i.e., dynamics 
governed by Pick's Law = D(j)"). Our experiment is 
also the first to observe persistence in one dimension (1- 
D). We employed a novel NMR technique to create an 
initial "quasi-random" spatial variation in the nuclear 
spin magnetization of a sample of laser-polarized -'^^^Xe 
gas. Subsequent 1-D NMR imaging allowed us to moni- 
tor the temporal evolution of the ensemble. We observed 
persistence by noting mean magnetization sign changes 
at fixed locations of constant size (i.e., imaging pixels) 
as a function of time. Using a simple theory (the "in- 
dependent interval approximation") and numerical sim- 
ulations, both Derrida et al. and Majumdar et al. ||] 
independently found that 9 w 0.121 for 1-D diffusion. 
Newman and Toroczkai found 9 « 0.125 in 1-D us- 
ing an analytic expression for the diffusion persistence 
exponent. Our measurements are consistent with these 
calculations. 

TABLE L A sample of reported persistence exponents. All 
values except those indicated are derived from numerical sim- 
ulations; (*) denotes exact analytical results, {}) experimental 
measurements, and (*) the result reported here. 

Dim. Diffusion Ising q-Potts 



1 0.12, 0.118* 3/8*, 0.35 

2 0.19 0.22, O.igt 

3 0.24 0.26 



(^) 
0.86, 0.88* (large q) 
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Recently, laser-polarized noble gas NMR has found 
wide application in both the physical and biomedical sci- 
ences. Examples include fundamental symmetry tests 
p2t , probing the structure of porous media |^^, and 
imaging of the lung gas space |Q . These varied investi- 
gations, as well as the experiment reported here, exploit 
special features of laser-polarized noble gas: the large 
nuclear spin polarization ('^ 10%) that can be achieved 
with optical pumping techniques; the long-lived nuclear 
spin polarization of the spin-1/2 noble gases ^^^Xe and 
•^He; and rapid gas-phase diffusion. 

We performed laser-polarization of xenon gas using 
spin-exchange optical pumping We filled a coated 

cylindrical glass cell (~ 9 cm long, 2 cm I.D.) with 
approximately 3 bar of xenon gas isotopically enriched 
to 90% ^^^Xe, 0.5 bar of N2 gas, and a small amount of 
Rb metal. We heated the sealed cell to ~ 100°C to cre- 
ate a significant Rb vapor. Optical pumping on the Rb 
Dl line was achieved with 15 W of circularly-polarized 
795 nm light (FWHM ~ 3 nm) from a fiber-coupled laser 
diode array. After 20 minutes the ^^^Xe gas was routinely 
nuclear spin-polarized to 1% by spin-exchange collisions 
with the Rb vapor. We next cooled the cell to room 
temperature in a water bath and placed the cell inside 
a homemade RF solenoid coil (2.5 cm diameter, 15 cm 
long, Q ^ 900) centered in a 4.7 T horizontal bore mag- 
net (GE Omega/CSI spectrometer /imager) with ^^^Xe 
Larmor frequency = 55.345 MHz. To allow the gas tem- 
perature to reach equilibrium, we left the cell in place 
for 20 minutes before starting the persistence measure- 
ments. At equilibrium under these conditions, the ^^^Xe 
polarization decay time constant (Ti ) was in excess of 3 
hours, with a ^^^Xe diffusion coefficient of 0.0198 cm^/s 
pTf . (Note that changes in the sample gas pressure, and 
hence the diffusion coefficient, simply cause a rescaling of 
the time variable and do not affect the persistence power 
lawp(t)^i-« 1.) 

The NMR pulse sequence we used to observe per- 
sistence in laser-polarized ^^^Xe gas diffusion is shown 
schematically in Fig. |l|. The first portion of the pulse 
sequence encodes a 1-D "quasi-random" pattern on the 
transverse magnetization of the laser-polarized ^^^Xe gas 
sample by using m pairs of variable-strength RF and 
magnetic-field-gradient pulses, repeated in rapid suc- 
cession (see Fig. 0). Each pair of RF and gradient 
pulses adds different spatial Fourier components to the 
1-D transverse magnetization pattern, with wavenum- 
bers given by the gradient pulse area and Fourier compo- 
nent amplitudes set by the RF pulse area. Next, a 7r/2 
RF pulse "stores" this quasi-random 1-D magnetization 
distribution along the longitudinal (z) direction while a 
subsequent strong (crusher) gradient pulse dephases any 
remaining transverse magnetization. The quasi-random 
longitudinal magnetization distribution then evolves with 
time due to diffusion and is monitored by a series of 1-D 
FLASH (Fast Low Angle SHot) NMR images MM (see 
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FIG. 1. NMR pulse sequence used to encode a 1-D 
"quasi-random" pattern on the average magnetization of 
laser-polarized ^^^Xe gas. Temporal evolution of the mag- 
netization pattern is monitored with n repetitions of a 1-D 
FLASH imaging routine j28|. For example, with m = 8 en- 
coding RF pulse/gradient pairs, the encoding pulse angles at 
= [30°, 35°, 37°, 41°, 45°, 50°, 63.5°, and 90°] while the 
encoding gradient amplitudes gi were chosen randomly. The 
imaging pulse angle Pj was fixed at 8° and the diffusion times 
Tj were varied from 2.4 ms up to ~ 2 s. The encoding gra- 
dients and the transverse magnetization dephasing "crusher" 
gradient were pulsed for 1 and 20 ms, respectively. Imaging 
gradients were applied for a total of 7.56 ms. The maximum 
gradient available was 6.7 G/cm. 



Fig. 0). 

The initial pattern of longitudinal ^^^Xe magnetiza- 
tion is quasi-random in that there must be a minimum 
length scale to the induced variations in the ^^^Xe mag- 
netization, i.e., a maximum wavenumber in the pattern, 
for there to be sufficient NMR signal for useful imaging. 
(This minimum length scale was « 0.6 mm in our exper- 
iment.) Nevertheless, at longer length scales the induced 
pattern must be random enough that persistence behav- 
ior can be expected. Ideally, ((/)(x, O)0(a;', 0)) = 5{x — x')] 
however, calculations indicate that it is sufficient for 
the initial condition correlator to decrease faster than 
\x — x'\~^ for 1-D persistence We found that six 
to eight RF/gradient pulse pairs (m = 6-8) were opti- 
mal for the desired quasi-random 1-D patterning of the 
^^^Xe magnetization, m < 6 resulted in a pattern that 
was not sufficiently random, while m > 8 significantly 
reduced the signal-to-noise ratio (SNR) of the NMR im- 
ages. The requirement of m > 6 is supported by numeri- 
cal calculations in which we modeled the NMR encoding 
sequence and simulated the subsequent gas diffusion us- 
ing a finite difference first-order forward Euler scheme 
[Q,^: we found persistence behavior (i.e., p{t) ~ t~^) 
only when m > 6. The requirement of to < 8 was set by 
the available NMR signal (i.e., the finite ^^^Xe magneti- 
zation), the necessity of rapid data acquisition to avoid 
excessive diffusion during the imaging sequence itself, the 
limitation of approximately 27rx (0.6 mm)~^ for the max- 
imum wavenumber, and the maximum imaging gradient 
strength available. 

For NMR imaging, we used a field of view (FOV) of 
31.5 cm with 0.6 mm resolution, which thus divided the 
9 cm cell into about 150 imaging pixels, each correspond- 
ing to a discernible spatial region of fixed size and loca- 
tion. We typically employed 8° excitation RF pulse an- 
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FIG. 3. A log-log plot of p{t), the fraction of regions (pix- 
FIG. 2. (a) Typical quasi-random initial pattern of '^'-"Xe els) with ^^^Xe spin magnetization that had not changed sign 



magnetization resulting from 8 encoding RF pulse/gradient 
pairs, (b) 32 images of the evolving magnetization pattern 
were acquired at logarithmically increasing times. White 
(black) indicates average positive (negative) magnetizations 
for each of the 150 fixed-location imaging pixels of width 
~ 0.6 mm. Gray indicates pixels whose magnetization has 
changed sign at least once. The persistence exponent is de- 
termined from the growing fraction of gray pixels as a function 
of time. 



gles and acquired 32 1-D images spaced logarithmically 
in time from ~ 3 ms to 5 s for a single experimental run. 
Longitudinal magnetization depletion due to imaging was 
highly uniform across the sample and did not affect the 
persistence measurement, since the relative magnetiza- 
tion amplitudes in neighboring imaging pixels was un- 
changed. An example of the images acquired in a typical 
run are shown in Fig. ^. For each pixel, we derived aver- 
age magnetizations (aligned or anti-aligned to the main 
magnetic field) from the phase information contained in 
the time-domain NMR image data, and spatial positions 
from the frequency information An experimental 

run thus provided a record of the ^^^Xe gas magnetiza- 
tion in each pixel as a function of time proceeding from 
the initial quasi-random pattern to the equilibrium con- 
dition of homogeneous (near-zero) polarization. 

To measure persistence, we noted the sign of the ^^^Xe 
magnetization in each fixed spatial region (i.e., in each 
1-D image pixel) and counted how many remained un- 
changed as a function of time. We equated the probabil- 
ity p{t) with the fraction of pixels that had not changed 
sign up to time t. We chose i — to coincide with the 
first image and assigned the time index for each image to 
be the start time of the imaging RF pulse. Images with 
SNR < 40 were excluded from the data to minimize un- 
certainty in pixel sign changes. We conducted about 30 
experiments with image SNR > 40, each with a unique 
set of randomly chosen encoding magnetic field gradi- 
ents {gi}. We observed that pixel sign changes occurred 
smoothly, so it was unlikely we missed sign changes 
with an error of more than one step in the imaging se- 



up to a time t, representing the sum of --^ 30 different experi- 
mental runs. The solid line is a weighted linear least-squares 
fit to the data for 0.1 s < t < 1 s, and yields 9 = 0.118±0.008. 
Error bars are derived from the number of pixels with ampli- 
tudes close to the image noise level and are shown when they 
exceed the plot symbol diameter. 



quence. We employed two averaging schemes to combine 
the results from different experimental runs. In the first 
method, we used a linear least-squares fit of log[p(t)] vs. 
log[i] for each run, resulting in a distribution of power 
law exponents with a weighted mean 9 = 0.119 ± 0.048. 
With our numerical simulations of quasi-random initial 
conditions, we found that this averaging scheme results 
in a gaussian distribution of exponents with a mean value 
6 w 0.12 in agreement with previous calculations for 1-D 
diffusion @|,|l| and our experimental results. In the 
second averaging scheme, we combined the data from 
all experimental runs; hence p{t) represented the frac- 
tion of total pixels from all experiments that had not 
changed sign up to time t. We found p{t) ~ with 9 
= 0.118 ± 0.008 for 0.1 s < t < 1 s. Figure | shows a 
log-log plot of p{t) vs. t when the data is averaged using 
this method. 

The observed deviations from power law behavior for 
t < 0.1 s and < ^ 1 s are explained by image res- 
olution and finite size effects, respectively. At short 
times persistence is not observed because '^^^Xe atoms 
have not yet diffused on average across a single 1-D im- 
age pixel (« 0.6 mm). The relevant diffusion time is 
(0.6 mm)^/(2_Djs:e) ~ 0.1 s for our typical experimental 
conditions. At long times, the coarsening of the ^^^Xe 
magnetization pattern by diffusion results in large "do- 
mains" of adjacent pixels with the same sign of the mag- 
netization. Our simulations indicate that persistence is 
not observed when there are few domains in the finite 
size sample because the number of magnetization bound- 
aries is greatly reduced; hence the rate of pixel sign- 
changing (i.e., the growth of the gray area in Fig. H) 
is curtailed. Both the short and long-time deviations 
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FIG. 4. The average size L of domains of adjacent pixels 
with uniform magnetization orientation, as a function of time 
t, derived from all experimental runs. For 0.1 s < t < 1 s, 
L ~ where ol = 0.45 ± 0.02 (solid line). The dotted line 
shows the expected L ^ t^^"^ behavior for an infinite system. 
The error in L is shown where it exceeds the plot symbol size. 
The finite size limit on L is evident in the four late-time points 
(A), which were taken from the only two runs with sufficient 
SNR at long times. 



are seen in Fig. ^ where the average length L of hke- 
signed magnetization domains from all experimental runs 
is plotted against time. For 0.1 s ^ t ^ Is, our data 
are in reasonable agreement with the expected power law 
L ^ t^/^ for diffusion For t > 1 s, we find L^l cm, 
implying there are typically less than 10 magnetization 
boundaries in the 9 cm long sample cell. 

In conclusion, we experimentally measured a persis- 
tence exponent 9 k, 0.12 for 1-D diffusion, consistent with 
analytical and numerical studies. We performed the mea- 
surement using a novel NMR scheme with laser-polarized 
^^^Xe gas which allowed us to both encode a spatially 
"quasi-random" magnetization pattern and monitor its 
evolution over several seconds. We also observed the ef- 
fect of finite sample size for long-time diffusion. In future 
work the experimental technique employed in this study 
may allow measurements of persistence in 2 and 3-D dif- 
fusion, in heterogeneous systems (e.g., porous media) in- 
fused with noble gas, and in 'patterns' p2[ . 
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